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Quantum graphity offers the intriguing notion that space emerges in the low energy states of the 
spatial degrees of freedom of a dynamical lattice. Here we investigate met ast able domain struc- 
tures which are likely to exist in the low energy phase of lattice evolution. Through an annealing 
process we explore the formation of metastable defects at domain boundaries and the effects of 
domain structures on the propagation of bosons. We show that these structures should have observ- 
able background independent consequences including scattering, double imaging, and gravitational 
lensing-like effects. 
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I. INTRODUCTION 

A single theory that reconciles quantum mechanics 
with general relativity would revolutionise our fundamen- 
tal understanding of reality. Such an endeavour in the 
search for a quantum theory of gravity has traditionally 
been approached either through the quantization of gen- 
eral relativity or extensions of quantum field theory. For 
example loop quantum gravity [1], spin foam models [2], 
causal dynamic triangulation (CDT) [3], and group field 
theory 0,[5| fall into the first category and collectively the 
string theories fall into the second [6]. Condensed mat- 
ter physics provides a third conceptual framework, where 
gravity is viewed as an emergent phenomenon. Motivated 
by techniques and concepts native to the study of many- 
body interactions in condensed matter physics, this per- 
spective has recently stimulated investigations into grav- 
ity analogues in condensed matter systems [p-[l2|, emer- 
gent graviton excitations in qubit model s fl3| , a nd quan- 
tum graphity (CDT and matrix models [l3:[lll may also 
be viewed from this perspective). 

Motivated by the removal of presumptions of the na- 
ture of spacetime, quantum graphity (QG) was proposed 
as a model in which ideas such as continuity, dimension- 
ality, and macro locality of the spacetime manifold are 
emergent phenomena. Previous works have studied how 
locality emerges in the model [16|, the role matter plays 
in the emergence of extended geometries [17], the entan- 
glement of matter with spatial degrees of freedom 1 181 . 
Ising mappings to study low temperature properties [19| , 
and the entrapment of matter in regions of high connec- 
tivity [20]. See Ref. [ITj for a review. 

In QG it is purposed that a low dimensional regular 
graph, representing fiat local space, emerged from an 
early universe, represented by a complete graph. As ev- 
ery vertex is connected to every other, there is no no- 
tion of subsystems and hence no notion of locality, which 
Konopka et at. [l^] interprets as representing a state with 
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no space. The evolution from the complete graph of the 
early QG universe to a lower energy state of the spatial 
degrees of freedom is reached by the destruction (and cre- 
ation) of edges. In the open (non-unitary) QG model [I3, 
it is assumed that the lattice is connected to an exter- 
nal heat bath through which edges are exchanged. Under 
particular parametrical constraints it was shown that the 
model can give rise to an hexagonal or honeycomb lattice 
as a stable local minimum [16]. This graph has the de- 
sirable properties that the spatial degrees of freedom are 
local and low dimensional. Here we explore metastable 
defects to this local minimum state and investigate the 
possible observable effects on the propagation of bosons. 

The study of defects is important, as in most real- 
istic or non-idealised systems, ranging from ferromag- 
net to field theoretic particle and cosmological models, 
metastable defects in some form exist. In many field 
theoretic early universe models for example, topologi- 
cal defects are an unavoidable causal process, and have 
been proposed to play an important role in cosmic evolu- 
tion [22|- In particular, cosmic strings have been pro- 
posed as the seeds for large-scale structure formation 
such as galaxies, offering a theoretical alternative to infla- 
tion [23|, |2^ . However observational data indicating that 
cosmic microwave background anisotropics significantly 
differ from the anisotropics that are predicted to be the 
result of topological defects, may relegate defects to a 
subsidiary role in cosmic structure formations [22]. In 
large extra dimensional models, topological defects pro- 
vide a mechanism by which 3+1 dimensional branes exist 
in a higher dimensional bulk |25l-[27|. 

Topological defects, of which the above are examples, 
occur as the result of the spontaneous choice of an order 
parameter value that breaks a symmetry of the system. 
Defects in QG on the other hand, of the type that are 
disruptions to the ordered crystal structure, are more 
akin to crystallographic defects. In this context, QG 
can be treated analogously to crystallographic models of 
condensed matter physics. Of course there are impor- 
tant differences. In conventional crystallographic models, 
the structures' internal inter connectivity is determined 
by the spatial arrangement of the atoms, which changes 



with the motion of the atoms. QG has an almost anti- 
thetical perspective: It is the interconnectivity of spatial 
points that determines spatial separation. The dynam- 
ics of the QG lattice is not viewed as the movement of 
these spatial points, but the change in the connections 
that relate these points. Moreover these interconnectivi- 
ties in QG are quantum degrees of freedom as they can 
be in a superposition of on or off states. Notwithstand- 
ing these differences, one can draw similarities between 
the evolution of the QG universe and the crystallization 
of solids. 

The ordering of unit cells into periodic structures (crys- 
tallization) is a well-known phenomenon. Although oc- 
curring in a wide variety of systems (e.g. freezing wa- 
ter, cooling magma), the process is always essentially 
the same: Random atomic distributions nucleate ordered 
structures below some critical temperature. And if this 
low temperature is maintained, these seeds will continue 
to grow and coalesce into larger crystal structures. The 
final atomic arrangement will be dependent on the initial 
configuration prior to cooling, rate of cooling, and ran- 
dom effects of thermal fluctuations. Typically a slower 
coohng rate will produce larger grains. Due to imperfec- 
tions and thermal fluctuations however, for most materi- 
als it is unlikely that the process will end in a structure 
without defects (although under very controlled environ- 
ments single crystals, such as monocrystalline silicon, can 
be manufactured). Instead the material will most prob- 
ably settle into a metastable (local minimum) state with 
many grains, resistant to thermal or mechanical pertur- 
bation. 

One may view the complete graph of the early uni- 
verse analogously to the high temperature diffused state 
of heated crystalline structures, and subsequently see in 
the QG model similar recrystallization qualities in the 
cooling phase. Specifically, fluctuations may seed the 
nucleation of space-like separated local energy minimum 
regions which are subgraphs of the ground state. The 
growth of these regions will give rise to domains (grains). 
As with the cooling of heated solids, it is not inevitable 
that the domains will coalesce to a global ground state. 
Analogous to conventional crystallization processes, the 
ripening (growth) of the different domains will see a com- 
petition for energetically favourable local configurations, 
resulting in a granular structure of space. 

Defects in the space manifold will affect particle prop- 
agation. On cosmological scales the classical treatment 
of the influence of gravity on matter has been sufficient in 
yielding accurate predictions: from the deflection of par- 
ticles near massive objects to gravitational redshifts. The 
fact that on large scales classical treatments accurately 
account for observation provides an important avenue to 
test QG. In particular, here we explore the possible obser- 
vational consequences of domain structures and defects 
in a classical lattice for the propagation of particles. 

In Sec. ini we briefly review and define the scope of 
the QG model to be considered. In Sec. IIIII we inves- 
tigate the formation of metastable domain structures. 




FIG. 1. A graph representation of the state vector 
1/12/13/14/23/24/34) = 1 101101) which has for 4 vertices and 
edges. 



Through an annealing process we show that defects that 
arise from an unstable interface of antiphasing domains 
are metastable. In Sec. IIVBI we examine the effects 
these structures have on boson propagation. We further 
explore other metastable domain structures that reveal 
interesting effects, leading to possible observable conse- 
quences. 



II. MODEL SETUP 

The Hilbert space of the spatial degrees of freedom 
in the QG lattice is the tensor product of the individ- 
ual state vector edge \lrs) ^ {|0), |l)}r,s- |0) r,s indicates 
no edge and |l)r,s an edge between vertices r and s. 
The set spanned by {|0), |1)} constitutes an orthonor- 
mal basis. We define creation and annihilation opera- 
tors on this space as b^^ = |l)(0|r,s and brs = |0)(l|r,s 
respectively. Fig. [1] provides an example of a graph 
representation of the state vector I/12/13/14/23/24/34) = 

IOl2|Ol3|Ol4|023|024|/)34) = |101101) in this SpaCC. 

Note only the interconnectivity of the graph is encoded 
in the state vector. The fact that we have embedded it 
in a two-dimensional plane and represented it as a square 
is only an illustrative choice. 

Matter degrees of freedom exist on the vertices of the 
graph. We consider the case of bosons. The Fock space 
is spanned by their individual vertex Fock states |n^). 
Acting on this space are the standard bosonic creation 
and annihilation operators, al and a^. The total Hilbert 
space consisting of the spatial and matter degrees of free- 
dom is spanned by {nr<s l^rs)(^nr' l^r')}- States of the 
system are represented as a tensor product of the spatial 
W = T^i^iYlrKsl^rs) and matter \(j)) = En^nllrl^r) 
degrees of freedom: \iIj) (g) \(j)). 

The edge number operator, which counts the number 
of edges at a vertex is defined as rrir = X^^ bl^brs and the 
particle number operator is defined in the standard way, 
rir = a^ar. Powers of the edge number operators are also 
defined as 
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Defined this way, rrirs gives the number of paths between 
vertices r and s of length L. 

In principle there are many forms the Hamiltonian H, 
which associates energies to the lattice, can take. How- 
ever the choice of H is greatly constrained by the QG 



proposition that locality is an emergent phenomenon. 
Konopka et al. ^16,] showed that a Hamiltonian consisting 
of valence term, i^vai, that energetically favors vertices 
with a parameterized number of edges and a loop term, 
H\oop^ that discriminately assigns energy to the number 
of closed loops, can give rise to a stable local minimum 
state that exhibits locality. In particular, the valence 
term is given by 



H. 



val 



gv 






^p{vo-mr) 



(2) 




b) 







c 












b' 


'""a 
4 




'd 


k 





where p is a real positive parameter that penalizes vertex 
valences (i.e. number of attached edges) that are away 
from vq^ and gy is assumed to be positive. The loop term 
is given by 
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(3) 
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where the coupling parameter gp is positive, e" 
^T=o ^^s sums all weighted pathways between ver- 
tex r and s; together with 5rs^ only closed pathways (or 
loops) are counted. The role of parameter R is to influ- 
ence the loop length L that corresponds to the peak of 
the effective coupling parameter, ^p^ = gpR^/L\. Specif- 
ically, ^p^ increases with L for small L, but quickly de- 
creases for large L; R determines this turnover point. It 
is of note that as the effective coupling rapidly decreases 
with large L, the contribution from i^ioop can be ap- 
proximated by truncating Eq. ([3j) at some maximal loop 
length; this is important as the loop counting process is 
computationally intensive. 

For the parameter choices vq = 3 and i? > 7.1 in the 
regime gy ^ ^p, J^vai + ^loop produces a local mini- 
mum honeycomb graph. Fig. [2^) shows a representa- 
tion of the honeycomb and Fig. [2)3) an isomorphically 
equivalent graph which we will call the brick representa- 
tion. Hys,i + i^ioop is by no means the only Hamiltonian 
that will result in low energy graphs with desirable traits 
such as locality and translation symmetry; other choices 
may yield different ground states also with such prop- 
erties. Without loss of generality, we use H^g,\ + i^ioop 
{vo = 3,i? > 7.1, ^v ^ ^p) as a toy example in a general 
framework that can be applied to other Hamiltonians, to 
gain insight into defects in the QG model. 

The evolution from the complete graph of the early 
QG universe to a lower energy state of the spatial degree 
of freedom is reached by the destruction and creation 
of edges. The lattice may be assumed to be in contact 
with an external heat bath through which edges are ex- 
changed [i^ or in an alternative unitary model, energy 
is conserved through the coupling of edges to matter de- 
grees of freedom [18]. Specifically in this unitary model, 
the destruction of an edge is accompanied by the creation 
of two bound quanta. The quanta hop according to 
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FIG. 2. Two isomorphically equivalent representations of 
a local minimum graph state for vo — ?>^ R > 7.1: a) 
Honeycomb representation [a — (0,0), 6 = ( — l,0),c = 
(0.5, 0.87), d = (2,0),e = (0.5,-0.87)] and b) brick repre- 
sentation [a = (0,0), 6 = (-1,0), c = (0,1), d = (1,0), e = 
(0, —1). a, 6, c, d, e are examples of the respective graphs' ver- 
tex labelling schemes. 




FIG. 3. A graphical representation of edge hopping term 
Hx — gx^q^gTnrsblqbqs. lu this cxamplc edge \lrq) is de- 
stroyed and edge \lqs) created. This local interaction is only 
allowed because there is a edge between vertex r and s i.e. 

TTT/rs I t-rs / — J- . 



where k, is proportional to the hopping frequency and 
rrirs restricts hopping to nearest neighbours. Note that 
in previous work [18| , t is used to denote this coupling pa- 
rameter. We use K instead, as t will be used to represent 
time. The bound quanta are then subsequently destroyed 
creating an edge to form a different lattice topology. 

Other possible dynamic interactions include edge ex- 
changes where the valence of vertices are conserved [1^, 
17] and edge hopping. An edge analog of i^hop, edge 
hopping is a local propagation of edges. We introduce it 
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where gx is a coupling parameter. Fig. [3] illustrates this 
interaction. In the context of the Hamma et al. [18,] 
model where every edge is assigned some constant en- 
ergy [/, every matter quantum assigned energy /i, and 
the edge-matter coupling constant is /c, Hx can be con- 
sidered as an effective dynamic lattice term in the limit 
/c <C |/7 — /i|. In this large detuning limit where the dif- 
ference in the lattice and matter energy scales are much 
greater than their coupling, the dynamics of the matter 
and lattice degrees of freedom are effectively decoupled. 
Here we investigate defects in the brick lattice. We 
study the large detuning or decoupled limit where matter 
dynamics (i^hop) are treated independently of the lattice 



dynamics (Hx) i.e. we ignore the effects of matter-lattice 
coupling. Furthermore, we restrict our investigation to a 
semi-classical model, in the sense that lattice dynamics 
are viewed classically through a heuristic process but the 
matter degree of freedom is studied quantum mechani- 
cally. 



III. METASTABLE DOMAIN BOUNDARY 
DEFECTS 

In the QG model, the early high energy high tempera- 
ture universe is represented with a complete graph which 
may be interpreted as a state without a notion of space. 
It evolves to a low energy low temperature state that 
exhibits translational symmetry and locality. The evo- 
lution process from the complete graph to this lower en- 
ergy state however is difficult to calculate as it amounts 
to a many-body problem for which the simulation time 
grows exponentially with the number of edges and ver- 
tices. Further compounding the situation is the problem 
of counting the number of closed loops. As the number of 
loops grows exponentially with the connectedness of the 
graph, this problem becomes increasingly difficult nearer 
to the early universe. 

Although a direct analysis of the evolution of the initial 
complete graph is impractical for graphs large enough to 
exhibit internal geometry, we can garner insight by com- 
parison with well known crystallographic systems. The 
random distribution of atoms, say in metal alloys, in 
the disordered high energy high temperature state means 
that the system is rotationally symmetric. Correspond- 
ingly for a complete graph with N edges, the system ex- 
hibits a discrete N -rotational symmetry [20|, which ap- 
proaches a continuous rotational symmetry as iV ^ oo. 

As the rotational symmetry of the condensed matter 
system which exists at high temperatures (energy) is bro- 
ken as the system cools, this also occurs in QG. Specif- 
ically, nucleation sees the ordering of atoms which grow 
into crystal structures that break rotational symmetry; 
and as the result of thermal fluctuations and other im- 
perfections, this crystallization process may occur inde- 
pendently in space-like separated regions. It is likely that 
a similar process occurs in the evolution of the complete 
QG graph to a lower energy state. It is however unclear 
whether quantum or thermal fluctuations or an interplay 
of the two will play the dominant role in the nucleation of 
low energy symmetry breaking in such crystalline struc- 
tures. 

In condensed matter systems, when there is enough en- 
ergy to overcome the activation barrier, the coalescence of 
locally ordered regions can allow the formation of larger 
crystal structures. However, even with a slow cooling 
rate, defects are likely to form in finite time. A par- 
ticular type of defect known as an antiphase boundary 
defect, occurs when the intersection of two domains are 
out of phase. To illustrate this, consider the crystalline 
structures of alloys which are interpenetrating lattices of 



b) 



FIG. 4. a) Two domains form independently. Where they 
meet is an unstable antiphase boundary defect, b) After 
the quenching process the boundary defect is frozen into a 
met ast able amorphous state. 



the different constituent atoms. For example alloys of 
two types of atoms, A and 5, of composition AB (e.g. 
CuZn) can form crystal grains of a simple cubic Bravais 
super lattice. This super lattice is composed of two inter- 
penetrating cubic sublattices of A atoms and B atoms. 
At grain intersections these sublattices may be out of 
step or phase, in which case an antiphase boundary de- 
fect occurs. 

Analogously domains of subgraphs of the ground state 
may also be out of phase at domain boundaries as de- 
picted in Fig. 3^). We will also consider examples of when 
domains differ in orientation in Sec. II V Bl and Sec. IIVCI 
Note that the metal alloy system involves the interleav- 
ing of two elements (e.g. Cu and Zn), whereas in the QG 
brick lattice the interleaving is purely a structural one, 
i.e. consisting of interleaving repeated structural layers 
as seen in Fig. |4^). For convenience we will refer to the 
brick graph as the ground state; whether it is actually 
the ground state or a local minimum does not matter 
for our purposes. The antiphase boundary represented 
by Fig. |4^) is unstable under lattice interaction. For ex- 
ample, edge hopping described by Hx will see an edge 
attached to a four edged vertex hop to a neighboring two 
edged vertex, producing a lower energy state of three 
edged vertices. In a full quantum mechanical treatment, 
complexity grows with the product of both the matter 
and lattice degrees of freedom. To simplify the problem, 
previous work [18] limited the manifold to four vertices, 
and the matter degrees of freedom were hardcore bosons. 
With such a small number of vertices, internal geometry 
is not a feature of the model, and so the lattice energy 
terms were simphfied to only being proportional to the 
number of edges. Here we are interested in the inter- 
nal geometry of the lattice and therefore must consider a 
manifold with many more vertices. The following section 
describes a heuristic approach to lattice evolution and 
formation of met ast able defects. 



A. Quenching 

We simulate the cooling of the lattice using the 
Metropolis algorithm [28]. The Metropolis algorithm [29] 
is an iterative algorithm where in each step the prob- 



ability P of evolving from state S with energy £^ to a 
neighboring state S' with energy E' is 
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where ks is the Boltzmann constant and T temperature. 
S' is randomly chosen from the neighborhood of states 
{S'} i.e. the set of states reachable by S. Specifically, 
motivated by the edge hopping term Eq. (|5]) , the dynam- 
ics of the lattice are mediated by the following heuristic: 
at each time step, h random number of edges are allowed 
to hop to their nearest neighbors, as depicted in Fig.[3l h 
corresponds to a measure of fluctuations. The states {S'} 
reachable by this interaction forms the neighborhood of 
states. 

As the probability distribution of states converges to 
the Boltzmann distribution for finite temperature on long 
enough time scales, Konopka et al. [17] used the Metropo- 
lis algorithm to study the ground state. Here we use the 
algorithm to study the stability of defects. 

We consider the scenario where domains which locally 
are ground states have formed, and that the tempera- 
ture has fallen below a critical value such that the sys- 
tem freezes into some local minimum. In particular we 
assume the extreme quenched case where the temper- 
ature is instantaneously reduced to zero. As the do- 
mains formed independently, their interface may be out 
of phase, for example as illustrated in Fig. |4^). This 
high energy configuration is unstable, and represents our 
starting state. To remove edge effects, we impose periodic 
boundary conditions, forming a torus. Along the domain 
boundary (major diameter) there are 100 vertices, and 
9 vertices running along the minor diameter. A small h 
will mean that the changes to the lattice are local pertur- 
bations, leading to the entrapment of met ast able defect 
states in local minima. Applying the Metropolis algo- 
rithm with T = 0, Fig. [IK) shows the decrease in lattice 
energy, averaged over 300 samples, at each t step of the 
quenching process. It shows that the average energy con- 
verges to 9.1 X l{fE/gY> and not the ground state energy 
8.8 X 10^^/^p. The parameter values are set to those 
that give rise to the brick graph as a local minimum, i.e. 
i;o = 3, r = 7.2, p = 1, gy/g^ = 10^. For practical 
computation times, we have limited the counted maxi- 
mum loop size, I/max = 10, in the above simulations. We 
have found no evidence that increasing Lmax significantly 
affects our results. 

There are many ways the lattice can evolve; moreover, 
under the extreme case of quenching the lattice will freeze 
into a state which is not the ground state. For example. 
Fig. [6^) shows a local part of a metastable configuration 
after quenching with h = 1. In this local area, the valence 
of vertex c is ^c = 2 and all other vertices have valence 
vo = 3. In the gy :^ gp regime, energetically nearby 
states are those where an edge has hopped to vertex c 
so that Vc = vq] all other neighboring states have more 
than one vertices with v ^ vo and hence are much higher 




t steps 



8.95 9.05 9.15 9.25 

E/gL,x10^ 



FIG. 5. Results of Metropolis algorithm applied to a toroid 
lattice (major diameter: 100 vertices, minor diameter: 9 ver- 
tices) with an unstable antiphase domain boundary [Fig. 3^)] 
at quenching temperature T = 0. Other parameters values: 
'yo = 3, r = 7.2, p = 1, gv/gp — 10^, Lmax = 10, iterations = 
300. a) Plot of the lattice energy averaged over 300 samples at 
each t step. The average energy converges to 9.1 x lO^^/^^p. 
This is higher than the ground state energy 8.8 x 10^ ^/^p 
(dotted line), meaning the lattice is more like to settle into a 
local minimum with defects rather than the ground state, b) 
The energy probability distribution of the lattice at t = 10^, 
when most of the samples have reached a metastable local 
minimum. As the lattice has been quenched, the distribution 
does not follow the Boltzmann distribution for finite T, but 
instead resembles the typical energy distribution of quenched 
metal alloys. 




FIG. 6. a) An extract of a metastable lattice at the domain 
boundary after quenching, b) Energetically nearby states ac- 
cessible from configuration a), by nearest neighbor edge hop- 
ping. The dotted lines represent hopping edges. Other neigh- 
borhood states (not shown) are of much higher energies as 
more than one vertices have v ^ vq. Calculation of if loop 
shows that the states in b) are higher energy states then con- 
figuration a). Therefore as configuration a) is energetically 
lower than its set of neighboring states, configuration a) forms 
a metastable state. 



energy states. There a four energetically nearby states, 
as shown in Fig. [6)3); however their ii'ioop is smaller than 
Fig. [6^), meaning that they have higher energies. There- 
fore the lattice freezes in configuration Fig. [6^), which 
is in a much higher energy state than the ground state 
due to the vertices which do not have valence vq. Note 
that if T > 0, then there is a finite probability of evolu- 
tion to higher energy states; in this case the lattice would 
eventually evolve to the ground state. 

The energy probability distribution of the local min- 
imum states reached by this quenching process is given 



in Fig. [SJd). As expected, it shows that the most prob- 
able state is centered around 9.1 x lO^E/gp. Treat- 
ing the problem as a statistical ensemble, we do not 
study the many possible individual metastable configu- 
rations. Nevertheless in general, the structures of these 
metastable states resemble two domains separated by an 
amorphous boundary defect. Fig.|4]3) and Fig. [6^) repre- 
sent extracts of two typical examples. In the local lattice 
region of configuration Fig. [4)3), all vertices have valence 
vo, therefore no further edge hopping is possible in this 
region, as doing so one would go to a higher energy state. 
In this local region. Fig. HJ)) is energetically separated 
from the ground state by the difference in loop contribu- 
tion, Hioop' As already discussed above. Fig. [6^1) is also 
metastable. Note that in the gv ^ gp regime, Fig.[6K) is 
a much higher energy state than Fig.|4)3). Inspecting the 
set of local minimum states of the numerical simulation, 
shows that there are a very low number of isomorphic 
graphs [30] . This indicates that the number of local min- 
ima in this energy landscape is large. As T = the dis- 
tribution is not the Boltzmann probability distribution of 
finite temperatures, but instead resembles the quenched 
energy distributions of, for example, amorphous metal- 
metaloid alloys [Slj. In the Boltzmann probability dis- 
tribution, lower energy states are more likely to form, 
with the ground states the most probable, which for the 
100 X 9 torus lattice corresponds to a ground state en- 
ergy of 8.8 X 10^ E/gp. In comparison, in the quenched 
probability distribution the most likely states are cen- 
tered around 9.1 x 10^ E/gp. In other words, under the 
prescribed conditions, stable local minimum defect states 
are more likely to form than the ground state. 



IV. EFFECTS OF DOMAIN STRUCTURES ON 
THE PROPAGATION OF BOSONS 

The presence of metastable defects affect the propaga- 
tion of particles. This is important as the observation 
of these effects could serve as a possible avenue through 
which to test QG. In this section we discuss the dispersion 
of bosons and simulate their propagation in the presence 
of defects. We study domain structures that give rise 
to scattering, refiection, refraction, and lensing-like be- 
haviour. 



A. Scattering 

In calculating the dispersion relation of bosons on the 
lattice we are only interested in the hopping term i^hop- 
We label the vertices of the ground state lattice with r^ = 
(xr^Vr) ^ N^. For a useful mapping between Euclidean 
space and the lattice, we associate a hop between near- 
est neighbors with one spatial u nit under the Euclidean 
distance function (i(ri,r2) = ^/{x2—Xl)^^^{y2—yl)^' 
This correspondence means that a neighborhood in Eu- 
clidean space corresponds to a neighborhood on the lat- 



tice. The brick representation as shown in Fig. [JJd) is a 
natural representation under this mapping and vertex la- 
belling scheme. By natural representation we specifically 
mean that embedding the lattice in Euclidean space, the 
vertex label corresponds to the Euclidean coordinates. 
As a comparison, a different labelling scheme, as indi- 
cated in Fig. [2^), sees the honeycomb vertex labels cor- 
respond to the coordinates of the embedding Euclidean 
space. In these examples, Euclidean space is the appro- 
priate embedding space as these lattices represent fiat 
space. However there are infinitely many possible embed- 
dings and each will have their corresponding dispersion 
relation. Contrasting such embeddings is not the aim of 
this paper, instead we examine scattering from defects 
in the brick representation. We however stress that scat- 
tering from defects is background independent, although 
the details of how they are represented are not, and so 
particle propagation is a useful tool for exploring poten- 
tial physically testable consequences of the QG model. In 
fact, detailed measurements of propagating fields may, in 
the future, allow us to reach some conclusions about how 
the microscopic structure of spacetime should be viewed. 
Forming primitive unit cells as shown in Fi g. [2b ), the 
matrix energy equation can be written as [32|, |33| 



Y,[Hrs]{<Pr} = E{cP,} 



(7) 



where [Hrs] is the (2 x 2) submatrix of the irreducible 
representation of i^hop that relates unit cell s with neigh- 
boring unit cell r ^ s\ [i^ss] is the intra-cell interaction, 
{^r} is a (2 X 1) vector denoting the wavefunction in unit 
cell r. Using Bloch's theorem, this equation can be solved 
with the ansatz 



K} = {0o}< 



ik D^ 



(8) 



where D^ invariantly translates the crystal from some 
arbitrary unit cell to unit cell r and k = (kx^ky) is a 
vector in the space reciprocal to D. 

Substituting Eq. ([8]) into Eq. ([7]) gives the energy eigen- 
value equation 



Y^lHr 



,ik-(D, 



'>{4>o} = E{M 



(9) 



which when solved gives the dispersion relation or band 
structure (see APPENDIX) 



[E{k.)/K]'^ = 1 + 4:Coskx{coskx + cos ky) . 



(10) 



This dispersion relation is notably anisotropic. However 
at low energies near k = the dispersion relation can be 
approximated to second-order as, ^(k) ^ ±/^(/c^ + /c^/3 — 
3) . This has the quadratic form of a conventional two- 
dimensional (anisotropic) free-space dispersion relation. 
The apparent anisotropy can be approximately elimi- 
nated by transforming the lattice (i.e. relabelling the ver- 
tices) to the isomorphic honeycomb representation. The 
dispersion relation in the honeycomb representation is 



£;hex ^ ±/^^l + 4cosV^A:a,/2(cosV3A:^/2 + cos3fe^/2) 
which in the low energy hmit is approximated as 
^hex(l^) ^ ±A^(3/c^/4 + 3^^/4-3), resembhng the disper- 
sion relation in conventional two-dimensional (2D) flat 
space. At higher energies however, higher order terms 
become significant implying that the dispersion relation 
in any representation is no longer quadratic, and the 
anisotropy can not simply be relabelled away. As we live 
in an isotropic universe (although there is some indica- 
tion that anisotropics exist on cosmological scales [34]), 
the honeycomb representation may be the appropriate 
representation in 2D. However to more clearly reveal the 
effects of the anisotropy of the lattice, we use the brick 
representation where the anisotropy of the dispersion re- 
lation is evident even at low energies. 

To see the effects of the antiphase boundary de- 
fects on the propagation of bosons, we simulate a sin- 
gle moded Gaussian pulse centered around (xo^yo) with 
width {ax.cTy), 



t=100/K 




FIG. 7. (Color online) Simulation of a localized boson as it 
scatters at a metastable amorphous domain boundary defect, 
a) The boson is initialized at xo = (100, 80), with cr = (10, 10) 
and k = (0, -2) so that it has v^/k = ±(0, 1). Most of the 
boson wavepacket propagates as vj" . Inset: Zoomed depiction 
of the amorphous defect line which runs from vertex (0, 100) 
to (100, 100). b) Snapshot of the system at t = 100/ k, showing 
that the amorphous defect has caused the boson to scatter; 
there is some probability that the boson is also localized on 
the domain boundary for a finite time. 



mO))=N^e- 



iix- 



■xof/2al]^-[iy-yof/2al]^ik^x^ikyy\'^\ 



(11) 

where A/" is a normalization factor. Evolution of the sys- 
tem is dictated by the Schrodinger equation, |V^(t)) = 
e*^'^°p^|'0(O)). It is important to note here that the evo- 
lution of the system is dependent on the lattice intercon- 
nectivity and not the representation. 

In Fig. [7^) a Gaussian pulse is initialized at xq = 
(100, 80) with (7 = (10, 10) and k = (0, -2). The gradi- 
ent of Eq. ([To]) gives us the group velocity (v^ = \/]^E). 
As we initialize the boson away from the domain bound- 
aries we can use the dispersion relation, Eq. (fTOj) . to esti- 
mate the group velocity, v^/a>: = ±(0, 1). Fig. (Tb) shows 
that most of the pulse being in the £^^(0, —2) state prop- 
agates as v+, and the metastable amorphous defect re- 
gion has a scattering effect on the boson. Although we 
have chosen to present the simulation in the brick rep- 
resentation, the manifestation of scattering from domain 
boundary defects are representation independent conse- 
quences of metast ability in QG. In other words, choosing 
a different representation changes only the details of the 
scattering, and hence scattered bosons may be used by 
an observer to detect the presence of a domain bound- 
aries. We also note that there is also some localization of 
the boson along the domain boundary for a finite time. 



B. Refraction and reflection 

The discrete ordered structure of crystals means that 
they are not rotationally invariant. In QG, crystalliza- 
tion of space-like separated regions into subgraphs of the 
ground state, having no preferred orientation, will give 
rise to dislocations when they meet. Fig. [8] illustrates 
such an example. The relative orientation of the domains 
can be represented by a global a e [0, tt) rotation of the 



a=0 



■■^1^ 



a=n/2 



domain boundary 



FIG. 8. Two domains formed independently with different 
orientation, labelled as domain- 7r/2 and domain-0 . At their 
intersection is a domain boundary defect. 



vertex labels. Within a domain there is short-range or- 
der, but this symmetry is broken over longer ranges. We 
investigate the effects of these types of domain structures 
on the propagation of bosons. 

As an introductory case we consider two domains that 
are oriented perpendicular to each other as represented in 
Fig. [HI This state is metastable under small fluctuations 
h as there is no nearby state of lower energy i.e. local 
edge hopping will not produce a lower energy state. We 
will identify domains by their orientation angle. 

The dispersion relation of the ground state with a gen- 
eralized rotation of the vertex labels is 

[£^(k, a)/K]'^ = 1 + 4 cos{kx sin a — ky cos a) 

X [cos{kx cos a -\- ky sin a) + cos{kx sin a — ky cos a)] . 

(12) 

Fig. [9] shows the isoenergy contour plot for a = and 



a^n/2 
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a=nf2 
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FIG. 9. (Color online) Isoenergy contours of the dispersion 
properties of domain-0 (left) and domain- 7r/2 (right). A mode 
with E/k = 2 (bolded contour lines) and ko = ( — 1.6, —0.5) in 
domain-0 has v^/a^: = (0.88, 0.83) as indicated by the arrow. 
Neglecting edge effects, phase matching will mean that this 
mode will couple to propagating mode \<.7^/2 = ( — 1-0, —0.5) in 
domain- 7r/2, which will have v^/^c = (1.6,0.3) (indicated by 
arrow) . 



a = 7r/2. 

Domain boundaries introduce edge effects, which 
means Eq. ([12]) is only valid away from the domain 
boundary. Nevertheless, as simulations will show, the 
dispersion relation can be used as a good indicator of 
the refractive and reflective properties across the domain 
boundary. 

The difference in the dispersion relation of domains 
of dissimilar rotation angles will mean that as the bo- 
son passes from one domain to another it will experience 
refraction. Incident and refracted angles are described 
relative to an interface normal. Using tan[^(k, a)] = 
Vg^y/vg^x aud Eq. ([12]) . the angle of propagation relative 
to the domain boundary normal is given by 



a=0 



a=n/2 




100 150 200 250 300 350 



400 



6>(k, a) = arctan 



CASBSa + {cb + 2ca)saCo 



SBCACa - {cb + '2CA)sASa 



, (13) 



FIG. 10. (Color online) Simulation of the propagation of 
a boson wavepacket over time. Different time instances, 
t, are superimposed; each instance is labeled with (/^t, M). 
For clearer presentation, populations are multiplied by factor 
M so that the peak value at each time instant is approxi- 
mately the same, a) The boson is initialized in domain-0 
with ko = (-1.6,-0.5). As the boson propagates through 
the domain boundary it undergoes refraction with angle of 
refraction \6r\ ~ 8.8^ . b) The boson is initialized in domain- 
7r/2 with k^/2 = (0.84, —1.1). As this mode can not couple 
to resonant modes in domain-0 the boson will be completely 
reflected. Inset: Zoomed depictions of the domain boundary. 



vation and phase matching can not simultaneously oc- 
cur and these field components will be reflected to form 
a second image of the object. Fig. [TOb) demonstrates 



where A = kyCa - k^s^ and B = k^Ca + kySa^ We have this with K/2 = (0.84, -1.1). The boson initialized at 



used the notation Cq = cosq and Sq = sing, q = A^B^ a. 

The refractive and reflective properties of the domain 
boundary are simulated in the propagation of a boson 
Gaussian wavepacket in a discrete space described by the 
geometry of Fig. [51 A boson initiated in domain-0 with 
ko^^o = (—1.6,-0.5) will have an incident angle given by 
Eq. ([T3]) as \6i\ = 43°. As the boson is initialized away 
from the domain boundary and hence edge effects, this 
propagation angle is well matched in the simulation as 
seen in domain-0 of Fig. [TOk). This boson will couple to 
a mode in the domain- 7r/2 with the same energy and ky 
component. Without edge effects, l<^7^/2 = (-1.0,-0.5), 
giving a refraction angle \0r\ = 8.8°. Edge effects at the 
domain boundary will cause deviation from this refrac- 
tion angle, however this angle can used as a guide as to 
how light will refract, as shown in Fig. [TOk). 

In this toy 2D universe an observer in domain- 7r/2 may 
see two images of an object. This is because for certain 
incident angles on the domain boundary energy conser- 



t = propagates hitting the domain boundary and gets 
reflected. 

We point out that although the calculated refraction 
and reflection angles are specific to the brick representa- 
tion and the domain structure setup in Fig. [8] the qual- 
itative effects of refraction and reflection, which are the 
result of the anisotropics of the lattice, will be present 
irrespective of the representation. In the hexagonal rep- 
resentation where the anisotropics are less severe than 
in the brick representation, refraction and reflection only 
become significant at high energies. Similar reflective ef- 
fects are also predicted with topological defects in cosmo- 
logical field theories [22]. For example the conical nature 
of space around straight sections of cosmic strings can 
give rise to double images of galaxies or quasars ^35| . If 
such observations are ever to be found, domain bound- 
aries may offer an alternative explanation. 

The observation of double imaging of astronomical ob- 
jects implies domain sizes at commensurate scales, how- 



ever it is not clear from our model what the average 
sizes of the domains should be. Although there is no 
precise method to predict the probable domain size in 
our model at present, it is likely related to the effective 
rate of cooling, as with conventional crystal properties. 
As fluctuations seed domains of random orientation in 
many different places, the rate of cooling will influence 
the ability of these smaller domains to coalesce into larger 
ones: Specifically a slower cooling rate would form larger 
domains than a faster rate. Many small domains with 
random orientations will form an amorphous manifold. 
With no preferred direction, on a scale much larger than 
the individual domains, this granular space will appear 
to exhibit rotational symmetry. On the other hand, at 
wavelength scales commensurate with the domain sizes 
the bosons will be scattered. In such a universe, dou- 
ble imaging of astronomical objects as a result of domain 
boundaries would not be observed. As our universe does 
appear to be isotropic and we do not see such scatter- 
ing effects, the length scales of the domains must either 
be very small (sub-microscopic) or very large (astronom- 
ical). 
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FIG. 11. (Color online) a) The intersection of four domains 
(bottom, left, top, right) with orientation a = 0,7r/2. b) 
Isoenergy contours at E/k = 2. The blue contour corre- 
sponds to domain-0 (bottom and top domains) and the green 
contour corresponds to domain- 7r/2 (left and right domains). 
The arrows represent group velocities. A mode originating 
in the bottom domain with k = (vr/G, — 7r/2) propagates with 
v^/k — (a/^/2, V3/2) (arrow 1). It couples to the propa- 
gating mode k = (0, — 7r/3) in the left domain (arrow 2). At 
the second domain boundary it will couple to the propagating 
mode k = (— 7r/2, — 7r/6) in the top domain (arrow 3). This 
trajectory will converge with its symmetric (about the y-axis) 
mode producing a lensing-like effect. 



C. Domain lensing 

Gravitational lensing is the bending of light as the re- 
sult of gravitational effects such that multiple or distorted 
images are formed of the source. Light is bent most near 
the center of the lens as the gravitational strength here 
is greatest. Because of this, gravitational lenses are char- 
acterized by a focal line (as opposed to a focal point of 
conventional lenses). Here we show how a granular struc- 
ture of space can give rise to similar lensing-like effects. 
This mechanism is distinct from the curved geodesies of 
general relativity. As an example we show that lensing- 
like effects can arise from the intersection of four (bot- 
tom, left, top, right) domains as illustrated in Fig. [TTk). 
We assume the structure to extend indefinitely. Because 
every vertex has valence vq = 3, the domain structure of 
Fig. [TTb ) is stable to small fluctuations {h = 1), as local 
edge hopping will not produce a lower energy state. 

To understand how this four domain geometry can pro- 
duce a lensing-like effect, we plot the isoenergy contours 
of the system at a particular frequency in Fig. [TTb). The 
blue isoenergy contour corresponds to domain-0 and the 
green to domain- 7r/2. The interface of the domains oc- 
cur at angles ±7r/4, requiring the phase matching con- 
dition kx 
terface and 



ky — 



fC™ 



k'^ — ky at the bottom-left domain in- 
ky = k'J. + ky at left-top domain in- 
terface with similar conditions at the right domain in- 
terfaces. For example a mode originating in the bot- 
tom domain with k = (— 7r/6, — 7r/2) propagates with 
v+//^ = (^/^/2, V^/2) (arrow 1). Phase matching at 
the domain boundary which is at 7r/4 relative to the x- 
axis, will mean it will couple to the propagating mode 
k = (0, — 7r/3) (arrow 2). At the second domain bound- 
ary which is at — 7r/4 relative to the x-axis, phase match- 



ing will mean that this mode will couple to the propagat- 
ing mode k = (— 7r/2, — 7r/6) (arrow 3). This example is 
symmetric about the x, y-axis. So together with its sym- 
metric mode reflected about the x-axis these trajectories 
will converge to focus the source. 

The domain lensing effect is simulated by propagat- 
ing a boson wavepacket in the lattice represented by 
Fig. [TTk). In Fig. [T2k). the boson is initialized as a 
Gaussian pulse of two modes, k = (±7r/8, — 7r/8), in 
the bottom domain. As it propagates through the do- 
main boundaries (indicated by the dotted lines) it un- 
dergoes refraction, following the predicted trajectories, 
to focus on the far side. The simulation of other modes 
[Fig. [T2b)x)] show that their foci do not converge to the 
same point, but instead form a focal line. It is important 
to point out that this behavior, where a point source is 
separated and then focused again on a contiguous line 
(as determined by lattice connectivity) is a property in- 
dependent of representation. In other words, this is a 
property of the interconnectivity of the domain struc- 
ture of Fig. [TTk ) and not it's brick representation. How 
we presented this in Fig. [12] is of course a result of our 
choice of representation. 

The presence of a focal line means that an observer in 
this 2D universe would detect a distorted image of the 
source. This effect can be compared to the distortion of 
images of astronomical objects due to gravitational lens- 
ing in our three dimensional world. It is important to 
note however that the effect in our model is purely non- 
relativistic. Furthermore, in GR gravitational lensing, 
light of all modes bends more the closer it is to the gravi- 
tational source; whereas in domain lensing (as illustrated 
in the preceding example) for a given mode the amount of 
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bending (refraction angle) is invariant with the distance 
from the optical axis. This distinction offers a possible 
avenue in which to test the model. 



V. CONCLUSION 

Quantum graphity is a background independent model 
that provides an alternative viewpoint on the notion and 
structure of space, based on condensed matter concepts 
but extended to a dynamic quantum lattice. Through an 
annealing process we explored metastable domain struc- 
tures and boundary defects. We investigated the propa- 
gation of bosons in these structures, revealing that they 
produce intriguing scattering, double imaging, and grav- 
itational lensing-like effects. Importantly this serves as a 
framework in which observable consequences of the QG 
model may allow it to be tested. 
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VII. APPENDIX 

Forming primitive unit cells as shown in Fig.[2l3) means 
that each unit cell will have four nearest neighbors. We 
label the displacement between a cell and its four neigh- 
boring cells as di = {-dx,dy),d2 = (dx.dy), ds = 
{dx^—dy)^d4 = {dx,—dy).The Hamiltonian in Eq. (J9j) is 
written as 



J2[Hr 



^ik-(D,-D,) 



oikd2 



ikd4 



1 



ikdi 



ikdc. 







(14) 



Solving for the eigenvalues of this matrix, the dispersion 
relation is 

[E{k.)/K]'^ = l-\-4:Cos{kxdy)[cos{kxdy)-\-cos{kydx)] . (15) 

In the brick representation (dx^dy) = (1,1) and in the 
honeycomb representation (dx^dy) = (3/2,>/3/2). 
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FIG. 12. (Color online) Simulation of boson Gaussian wavepackets with a) k = (ztyr/S, — 7r/8), b) k = (±7r/6, — 7r/6), c)k = 
(di7r/5, — 7r/5) at the four domain structure of Fig. 1111 The arrows indicate the group velocity (v^) in each domain. As the 
boson crosses the domain boundaries it undergoes refraction to converge on the far side. The foci points for the different modes 
of the boson wavepacket form a focal line. Note that there are also modes (the faint light blue pulses) associated with v^ which 
propagate in the opposite direction (i.e. downwards) and are reflected off the hard wall boundary. The notation (/^t, M) follows 
Fig.Iini 
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